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Mechanical properties are of central importance to materials sciences, in particular if they depend 
on external stimuli. Here we investigate the rheological response of amorphous solids, namely col¬ 
loidal glasses, to external forces. Using confocal microscopy and computer simulations, we establish 
a quantitative link between the macroscopic creep response and the microscopic single-particle dy¬ 
namics. We observe dynamical heterogeneities, namely regions of enhanced mobility, which remain 
localized in the creep regime, but grow for applied stresses leading to steady flow. These different 
behaviors are also reflected in the average particle dynamics, quantified by the mean squared dis¬ 
placement of the individual particles, and the fraction of active regions. Both microscopic quantities 
are found to be proportional to the macroscopic strain, despite the non-equilibrium and non-linear 
conditions during creep and the transient regime prior to steady flow. 


The properties of materials not only depend on their chemical composition, but also on the arrangement and 
dynamics of their constituents. It is thus crucial to understand the link between the macroscopic behaviour and 
the microscopic single-particle level. The relation between an applied mechanical force and microscopic processes is 
understood for crystalline, i.e. ordered, materials. Crystalline solids (like metals, ceramics or minerals) irreversibly 
deform when subjected to a load which is small enough to avoid fracture, although this response is very slow. This 
kind of response is called creep and originates from the presence of defects in the otherwise ordered arrangement of 
atoms. The diffusion of vacancies and dislocations is responsible for the observed plastic deformation [T]. The same 
relation and microscopic processes cannot occur in amorphous, i.e. disordered, materials. 

Nevertheless, in amorphous solid-like materials, a similar macroscopic creep response is observed under application 
of shear stresses below the yield stress CTy, i.e. below the transition from an elastic to a plastic response. The 
macroscopic creep response has been intensively studied in metallic, polymeric and colloidal glasses [M]. Several 
models mm successfully describe the time evolution of the strain measured during creep, namely its characteristic 
sub-linear time dependence. However, the relation of the creep response to the microscopic structure and dynamics 
has hardly been determined and is not well understood. Due to the disordered structure of amorphous solids the 
concept of defects is not applicable and a microscopic mechanism different from the one in crystalline solids must be 
responsible for creep. Thus, to make progress, microscopic observations on a single-particle level during creep tests 
are required. 

Combining experiments and simulations, we investigated colloidal glasses when constant stress is suddenly applied, 
i.e. during creep tests. We reveal a quantitative link between the macroscopic rheological response and the microscopic 
dynamics. This is possible due to recent developments in simultaneously performing rheology and confocal microscopy 
[igiiT]- During creep flow near the yielding threshold, we observe that very few particles undergo large non-affine 
displacements which leads to pronounced, but spatially localized, dynamical heterogeneities and sub-diffusive dynam¬ 
ics. In contrast, for stresses beyond the yield stress, transient super-diffusive dynamics mark the onset of steady 
flow. At the same time, growing domains of enhanced dynamic activity are present, with their number correlating 
with the macroscopic strain. This is reflected in a correlation between the macroscopic strain and the single particle 
displacements. In addition to the steady-state flow regime, this correlation also holds in the creep and transient 
states, specially for stresses near the yield stress. Hence, we can quantitatively relate the macroscopic rheological 
response of soft glasses to the average and heterogeneous microscopic dynamics which are spatially localized during 
creep but span the entire system at large stresses that lead to flow. The different microscopic behavior thus reflects the 
different macroscopic response during creep and flow, respectively. This extends previous observations to non-linear 
and non-equilibrium situations. Furthermore, as we observe the same behavior for different systems, realized in the 
experiments and simulations, this appears to be a general feature of glasses. 
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RESULTS AND DISCUSSION 

In our experiments and simulations we investigated two different model colloidal glasses. In the experiments, the 
glass is a binary mixture of sterically stabilised PMMA spheres with a size ratio of 5, dispersed in a density and 
refractive index matching solvent, with total volume fraction (j) = 0.61 and a relative volume fraction of small spheres 
xs = 4’s/'P = O-l- III this binary glass, the motion of the large particles is arrested via caging by neighbouring 
large particles p^l - l20] . In our molecular dynamics simulations, the glass is formed by a binary Yukawa fluid of 
large and small spheres with size ratio 1.2 and a relative number fraction of small spheres of 0.5, large enough to 
prevent crystallization. This system is quenched to T = 0.10, i.e well below the mode-coupling critical temperature 
of the system, Tc = 0.14. All times are normalized; in the experiments by the short-time one-dimensional Brownian 
diffusion time of the large spheres, = Sirrjd^/SkBT « 3.76 s, where du is the diameter of the large spheres, r] 
is the viscosity and feeT' the thermal energy, and in the simulations by the time unit Tq™ = (rndg/ess)^^^: where 
TO, ds and ess are units of mass, length and energy, respectively, with ds the diameter of the small spheres and ess the 
energy-scale corresponding to the interaction between small particles. The colloidal glasses investigated in experiments 
and simulations hence involve different interactions and different mixing and size ratios of their components. Using 
these different model systems allows us to explore the general features in the response of glasses to externally applied 
stresses. 

Macroscopic Strain is Quantitatively Related to Single-Particle Mean Squared Displacement. We 

performed a step to an applied constant stress (ct = const) on an initially quiescent glass. In the experiments, the 
stress was applied using a commercial stress-controlled rheometer, while in the simulations one wall was pulled by a 
constant force Fq. We monitored the macroscopic response via the time evolution of strain This situation is in 
contrast to the case of imposing a constant shear-rate (7 = const) pTll27) . where the bulk stress a{t) is monitored. 
Unlike for an applied shear rate 7 , when constant stress a is applied there is no timescale imposed and flow regimes 
below yielding can be investigated. The choice of control parameter, i.e. constant a or constant 7 , hence determines 
the intermediate flow states via which a glass evolves from the quiescent state to steady flow m- In the following, 
we exploit these possibilities and link the increasing macroscopic strain to the evolution of local particle motions, 
using stress as the external variable and including stresses below the yield stress Uy. In experiments, we estimated 
the yield stress of the glass, Uy « 0.010 Pa, from the stress at the crossing point of the storage and loss moduli in 
large amplitude oscillatory shear measurements at 1 rad/s. In simulations, at T = 0.10 the yield stress tJy = 0.072 (in 
simulation units) was estimated by strain-rate controlled simulations |28j . 

If the applied stress cr « Uy, a characteristic creep response is observed with the strain increasing sub-linearly 
with time within the experimental window, 7 ~ with a « 0.5 (Fig. broken line). Furthermore, for a = 0.9Cy 
(Fig.[T>, broken line), a smaller effective exponent is found, in agreement with previous results [51 [T211^51150] . Hence, 
for cr < cTy, the deformation occurs extremely slowly and the system is unable to reach a steady state within the 
observation time. This is reflected in the particle dynamics in vorticity (neutral) direction, namely the mean squared 
displacement (MSD) Ay^(t) (Fig. [l:,d, broken lines). In experiments and simulations, at short times the increase of 
the MSDs is limited, consistent with caging, while at longer times a sub-diffusive regime is observed; Ay^ ^ with 
6 < 1. We find 6 « a within the explored time window. The MSDs show little change with the waiting times after 
the beginning of the stress application (Fig. [^,d, broken lines). The observed macroscopic creep response and the 
absence of steady-state flow is thus connected to the particles’ inability to diffuse. 

In contrast, if cr ^ (Ty, the strain response shows a rapid transition to a steady flow regime, which corresponds 
to 7 ~ t, i.e. 7 is constant (Fig. [^,b, solid lines). The MSDs again display caging at intermediate times (Fig. [^,d, 
solid line), while at long times diffusion occurs [T7lI3TH33] . The slightly lower MSD plateau observed in experiments 
is due to cage constriction |18j and is also observed in Brownian dynamics simulations miiz], but not in molecular 
dynamics simulations where the microscopic dynamics is Newtonian. In between caging and long time diffusion, a 
transient super-diffusive regime is observed. This coincides with the transition of the rheological response from the 
initial elastic regime to the flow regime. Note that in the experiments, the initial superlinear increase in strain is a 
known effect of the rheometer’s inertia [M]. With increasing waiting time (Fig. [^,d), super-diffusion occurs at 
increasingly earlier times and for increasingly shorter time intervals, until it almost disappears in the steady state. 
Thus, the onset of flow is characterized by transient super-diffusion and, subsequently in the steady state, by diffusion. 
This indicates that the different regimes in the macroscopic strain response 7 (t) are reflected in different features of 
the single-particle dynamics, here characterized by the MSD Ay^(t). 

We now quantitatively investigate the relation between the macroscopic strain ■j(t) and the microscopic MSDs 
Aj/^(t). In the case of steady flow ^{t) = qf and, since then the particles diffuse, Ay‘^{t) ~ D{a)t, which implies that 
Ay‘^{t) ^ [D{a)/'y]j{t) = C{a)j{t). Previous experiments and simulations under constant applied shear rate have 
found D ^ ® [mill], which implies C'(cr) ^7 since stress and strain control are equivalent in steady flow. In 

our case, in the asymptotic diffusive regime (corresponding to 7 ^ 10, Fig. ,b) we observe an approximate linear 
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FIG. 1. Comparison of (left) experimental and (right) simulation results, (top) Time-dependence of the strain 7 ( 1 ) for applied 
stresses cr as indicated, relative to the yield stress ay. (bottom) Mean squared displacement in the vorticity direction, 
(indicated by same colors and line styles), immediately after stress application, i.e. for waiting time tw = 0 , and larger (as 
indicated) until the steady-state is reached, i.e. tw —>■ 00 (symbols). For the smaller applied stress, is divided by a factor 
3 for clarity, both in experiments and simulations. 


relation Ay^(t) ^ 7 ( 1 ) for a large range of a (Figj^, Sec. 1 in Supplementary Information). The slight shifts between 
the curves for different a occur due to the expected behaviour of C'(cr) (Fig. [^). If is rescaled by C{a), the data 
fall onto a single line of slope I (Fig. %)■ 

Although our argument for the relation Ay^(t) ^ 7 ( 1 ) is based on the assumption of steady flow, the relation sur¬ 
prisingly also holds in non-steady states for 7 < 10 , which corresponds to creep (for a ^ dy) or the transient regime 
before steady flow (for cr > cry). In contrast, Ay'^{t) ~ 7 ( 1 ) does not hold for large stresses cr and small strains 7 (or 
short times t). In both, experiments and simulations, systematic deviations are seen to occur with increasing stress. 
The deviations occur due to a time lag between the particles’ motion beyond their cages and the onset of macroscopic 
deformation (Fig. SM-1, Supplementary Information). Moreover, at very short times, i.e. in the initial elastic regime 
(Fig.[I]), the proportionality is also not observed. This suggests that the observed correlation is a consequence of 
the plasticity that develops after the initial elastic regime. Our observations mark a clear difference between the 
yielding response under applied constant stress (investigated here) and applied constant shear-rate (investigated in 
[2^ I26| ). In the latter, Ay^(t) ~ ^{t) cannot hold in the transient regime, where Ay^(t) increases superlinearly with 
t while 7 ( 1 ) increases linearly. Note that this connection between nonlinear strain and the single-particle dynamics is 
an implicit assumption in a recent theoretical approach based on a nonlinear Langevin equation |35H37] . Our data 
indicates to what extent such a connection is valid. 


Displacement Distributions Indicate Small Fraction of Mobile Particles. In addition to the characteri¬ 
sation of the particle displacements via the MSD, i.e. a mean value, we have also investigated the distribution of the 
displacements, namely the self part of the van Hove function p{Ay). For tr « cry (Fig. left), at all times the van 
Hove functions exhibit a nearly Gaussian shape for small Ay, which corresponds to localised particles, and moderate 
exponential tails which correspond to large displacements of a small fraction of particles. The non-Gaussian tails 
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FIG. 2. (a) Mean squared displacement in vorticity direction, Aj/^, as a function of strain 7 for different values of the applied 

stress (yjoy obtained in experiments and simulations (Sec. 1 of Supplementary Information). The experimental Ai/^ values 
are multiplied by a constant factor in order to match the simulation data, (b) Ratio of diffusion coefficient to shear rate, 
C{a) = D{a)/^, obtained from fits (Sec. 1 of Supplementary Information), as a function of the shear-rate in the steady state, 
7 . The dashed line indicates a power-law C{a) ~ 7 “°'^. (c) Scaling plot of Ay'^/C{a) as a function of 7 , with the dashed line 
indicating a slope of 1 . 


only slightly change with increasing time. This indicates that shear-induced delocalisation is a very slow process. In 
particular, large displacements at the shortest time of the measurement = 0.0037 hardly occur and therefore 

macroscopic how is delayed. 

For a > cTy (Fig. right) shear leads to a larger deviation from a Gaussian distribution with a signihcant number 
of large displacements. The deviation from Gaussian behavior was quantihed by the time dependence of the integral 
/R,(t) of the residuals of the Gaussian Hts to p(Ay), which, for each hxed time, was normalised to the integral of the 
distribution (Fig. inset). A non monotonic trend of /R(t) is observed, with a maximum value during the interme¬ 
diate superdihusive regime. At later times Init) continuously decreases and eventually vanishes when diffusion sets 
in and a Gaussian distribution of displacements is recovered. 


Evolution of Dynamical Activity Follows Macroscopic Strain. The tails in the van Hove function p{Ay) 
reveal the existence of a small number of very mobile particles during the transient regime. We quantify the time 
evolution of the fraction of these mobile particles by the ratio <i)^^(t) = /a(t) with the integral of p{Ay) for 
displacements Ay/d > SAymin, with d = ds and dh in simulations and experiments, respectively (Fig. solid 
symbolsb The value Aj/min = Ay"^(to)/d? is the localization length estimated from the MSDs at the shortest time 
to (Fig. [^,d). In simulations the time-dependence of $^^(t) closely follows that of the instantaneous strain 7 ( 1 ), 
up to 4>J^(t) = I (Fig. lines). In the experiments, similar results for are observed (Fig. [^) except that, 

in contrast to the simulations, ■jit) is not the instantaneous strain but a time average, leading to a small deviation 
between $l[^(t) and 7 ( 1 ). The macroscopic strain is therefore not only proportional to the average mean squared 
displacement (Fig.[^ but also the fraction of mobile particles: this indicates that the mobile, dynamically active parti¬ 
cles contribute most signihcantly to the mean squared displacement. This is true both below and above the yield stress. 


Spatial Distribution of Dynamical Activity is Heterogeneous. We introduce spatial coarse-graining in order 
to reduce noise. We divide the held of view into 10 x 10 square boxes, each with size (2.8 For each particle i, 

the displacement in the vorticity direction, Ayi{t) = yi{t) — yi(to), was determined. The average particle mobility in 
box Im, with l,m = 1... 10, was calculated according to 


= {Ayi{t))i^ 


( 1 ) 







5 


EXPERIMENTS 



SIMULATIONS 



Ay / ds Ay / ds 

FIG. 3. Van Hove self-correlation functions, i.e. distributions of displacements p{Ay), determined (top) by experiments for 
(a) a stress <j ~ ay, a waiting time tw = 0 and times = 0.0037, 0.54, 1.05, 2.12 and 5.31 (left to right) and (b) 

a ~ 5ay, tw = 0 and same times, except the longest time here is t/rQ^^ = 2.65, (bottom) by simulations for (c) a « l.luy, 
fw = 0 and t/TQ^ = 20,3.7 x 10^, 11 x 10^,56 x 10® and 110 x 10® (left to right) and (d) a = 1.53(Ty, tw = 0 and = 

20,1.01 X 10®, 3.7 X 10®, 7.4 x 10® and 18.6 x 10® (left to right). Dashed lines represent Gaussian hts to p{Ay) for small Ay. 
Insets: normalized integral 7 r of the residuals of the Gaussian fits in the main plots, as a function of time t/ro- 


where (... denotes an average over all the particles which were in the box Im at t = Iq. A box Im is defined active 
at time t if pim > SAj/min, following the same criterion used to distinguish largest displacements of single particles in 
the van Hove functions (Fig. ^,d). The fraction of active boxes, $^(0 = Na{t)/Ntot, with the number of active 
boxes and Ntot the total number of boxes. With time d)^(t) grows as the fraction of the single mobile particles $a^(t) 
(Fig. symbols). Thus, the time-dependence of d)^(f) is also proportional to 7 (t). A similar connection between 
the number of active regions and strain growth was experimentally observed in the creep flow of frictional granular 
particles [38] . 

To investigate the existence of heterogeneity in the dynamical activity, we consider the spatial distribution of active 
boxes. For cr « try, the distribution of local mobilities within the velocity-vorticity plane does not indicate any 
prominent features (Fig. top). At any specific time, there are some active boxes with larger mobilities, but the 
locations of the boxes with the largest mobilities vary randomly with time. For cr « ScTy, similar mobilities occur at 
short times, when the localisation plateau in the MSD is observed (Fig. [5}i,b, bottom). In contrast, at t > O.Stq’^^, 
roughly coincident with the onset of super-diffusion in the MSDs determined for ty, = 0 (Fig. a region with 
enhanced mobilities emerges (Fig. [^,d, bottom), expands with time (Fig. [^, bottom) and spans almost the whole 
field of view once the system flows (Fig.[^, bottom). Hence, the onset of flow (Fig. [I^,b) coincides with the appearance 
of a region of higher local mobility (Fig. and super-diffusive dynamics (Fig. |^,d). Furthermore, it leads to the 
pronounced non-Gaussian tails in the van Hove correlation function at intermediate times (Fig. which disappear 
once steady flow has developed and the dynamics again becomes more homogeneous (Fig.[^ ins^. 

The enhanced local mobilities do not result from sudden large displacements, but occur through the accumulation 
of only slightly above-average local, non-affine particle displacements. This has been confirmed by calculating the 
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FIG. 4. Fraction of active particles (solid symbols) and active boxes (open symbols) determined (a) by simulations 

at stresses o joy = 0.9 (•, o), 1.1 (a, A), 1.4 (0), 1.53 (>) and (b) by experiments at a jay ~ 1.0 (■,□) and 5.0 (A, A). Lines 
of the same colour represent the strain 7 for the corresponding applied stresses, where the instantaneous strain is shown in the 
case of the simulations. 



(a) (b) (c) (d) (e) (f) 


FIG. 5. Maps of average particle mobilities iiimit) within boxes Im (Eq[^ for (top) stress ct ~ o-y and (bottom) a ~ 5ay and 
times = 0.025,0.20,0.27,0.43,0.53,0.80 (a-f, indicated in Fig. by dashed lines) as observed in experiments. The box 

size is ( 2.8 di)^. 


instantaneous mobilities from O.lSr™’’ to 0.46 tq’'^ and 0.28Tg’'^ to 0.56 t™p, i.e. for 10 sampling times, instead of 
starting from the shortest measurement time (as in Fig. [5). No large instantaneous mobilities and no significant 
difference to cr « CTy are observed (Sec. 2 in Supplementary Information). Similar results are obtained in our simula¬ 
tions. The occurrence of correlated plastic events [ssmo] and avalanche-like behavior mma have been proposed as 
mechanisms driving the onset of flow. Such cooperative events might be connected to the correlated local mobilities 
and their spreading observed in our study. The observed intermittency in the displacements might also be related to 
stick-slip motion [15] . 

The larger area in the velocity-vorticity plane monitored in the experiments allows us to quantitatively investigate 
the spatial growth of active regions. If the box Im is active or inactive, nim is defined as I or 0, respectively. 
Based on this definition, we calculate the spatial correlation of active boxes, that is the box-box correlation function, 
G(r) = {nimnv^>) with = {l-l'f + {m-m'f (Fig.[^). The brackets (...) indicate an average over the individual 
boxes. The characteristic length ^ of the spatial correlation G{r) was determined by fitting a stretched exponential 
function /(r) = 7lexp[—(r/^)^] to G(r). The correlation length ^{t) increases from an initial value ^ « 5c?l at 
t = O.IOtq’^^ to ^ « SOdp at t = 0.92x0^^^, with ^{t) ~ (Fig. [^). For a ~ Uy the correlation length ^(t) instead 
does not grow and stays approximately constant for all times t (data not shown). 
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FIG. 6. (a) Box-box correlation functions G(r) for stress cr/o-y « 5 and time — 0.1, 0.20, 0.27, 0.43, 0.53, 0.66, 0.80 

and 0.93 (left to right) as observed in experiments. Lines represent stretched exponential fits, (b) Correlation length of active 
boxes, as a function of time for <j/ay « 5.0; the line indicates ~ 


CONCLUSIONS 


Using experiments and simulations, we demonstrated that under applied stress, the macroscopic deformation of 
glasses can be linked in a consistent way to the single particle displacements. In particular, the strain is approximately 
linearly related to the single-particle MSD even in the time-dependent non-linear response regime, including the creep 
and the transient regime preceding steady flow. Furthermore, the fraction of active particles in the van Hove function 
as well as the fraction of active regions, i.e. of groups of particles, is also proportional to the macroscopic strain. 
Heterogeneities in the location of these active particles are present both for applied stresses smaller and larger than 
the yield stress. The spatial distribution of regions with larger displacements determines the onset of flow. For 
applied stresses around the yield stress, i.e. during creep, localised regions of enhanced dynamical activity allow only 
for sub-diffusive dynamics. Increasing the stress beyond the yield stress, the active regions grow heterogeneously and 
super-diffusive transients emerge leading to particle diffusion with steady flow setting in. We observe qualitatively 
the same behavior for the different models studied in our experiments and simulations and thus expect that our 
observations represent generic features of glasses. 

Future work should focus on understanding how the external stress leads to the occurrence of locally enhanced 
mobilities, e.g. whether these are related to thermally activated local structural changes. Furthermore, the mechanisms 
that drive the spreading of the active regions within the plane as well as in the transverse direction need to be identified, 
thereby providing possible links to transient shear banding in the velocity-gradient plane [28l |44j . All these would 
help to develop a more complete scenario for the fluidisation of glassy systems under applied stress. Furthermore, it 
can open the route to the rational design of materials with desired response to applied stresses. 





METHODS 


Experiments. We investigated a mixture of sterically stabilized PMMA spheres of diameters dn = 1-76 /rm 
(fluorescently labeled) and ds = 0.35 /rm, dispersed in a cis-decalin/cyclohepty 1-bromide mixture which closely matches 
their density and refractive index. After addition of salt(tetrabutylammoniumchloride), this system presents hard- 
sphere like interactions [JSl IHj. The total volume fraction is 0 = 0.61 and the relative fraction of small spheres 
xs = 4‘s/4“ = 0.1. The formation of a glassy state in this mixture was demonstrated by using rheology and confocal 
microscopy measurements of the dynamics of large particles |18t42()j . The presence of small spheres, with their larger 
energy density, increases the yield stress of the system, thereby improving the quality of the rheological data while still 
allowing for the simultaneous observation of the large spheres with confocal microscopy nmn]. The rheological and 
confocal microscopy measurements reported in the manuscript were obtained using a combination of a commercial 
MCR-301 WSP stress-controlled rheometer (Anton-Paar) and a VT-Eye confocal unit (Visitech) mounted on a Nikon 
Ti-U inverted microscope, with a Nikon Plan Apo 60x oil immersion objective (NA = 1.40). We used a cone-plate 
geometry of diameter 50 mm, cone angle 1° and truncation gap 100 /rm. The bottom plate consists of a microscope 
coverslip which was coated with a mixture of PMMA particles of radius 0.885 /im and 0.174 ^m. The surface of the 
cone is sandblasted. The roughness of the geometries prevents wall-slip, as verified by imaging. A solvent trap was 
used to reduce solvent evaporation. Due to the fact that rheological measurements on colloidal glasses are affected by 
loading effects, shear history and aging, before each test a renjuvenation procedure was performed in order to obtain 
a reproducible initial state of the system. After loading, we performed a dynamic strain sweep to estimate the yield 
strain 7 y of the system from the crossing point of the strain-dependent storage, G", and loss, G", moduli. Oscillatory 
shear at 7 = 3 7 y was applied to induce flow and maintained until the time-dependent G' and G” reached a 
stationary state, typically after 200 s. Afterwards, oscillatory shear in the linear viscoelastic regime, 7 = 0.001, 
was applied until G' and G" became stationary, typically for t > 300s. The state characterised by the stationary 
values of G' and G” was the initial state, prepared before each creep measurement. The relative error on the strain 
determination during creep measurements is smaller than 1 %. 

Confocal microscopy images were acquired in a velocity-vorticity plane about 6 mm from the center of the geometries 
and 30 /im from the bottom plate. Images with 512x512 pixels, corresponding to 51 /tm x 51 /im, were acquired at 
a rate of 67 frames per second, which ensured accurate particle tracking even at the highest applied stresses (typical 
movies in Supplementary Information). By imaging the truncation gap of the cone, we verified that bending of 
the coverslip is negligible m- This is also indicated by the fact that, despite the applied stress, the particles in 
the imaging plane remain perfectly in focus (movies in Supplementary Information). The fact that we can image 
the truncation gap of the cone is also used to check that the bottom plate is perpendicular to the rotation axis of 
the cone. Particle coordinates and trajectories were extracted from the images using standard routines [49) . Mean 
squared displacements from four independent measurements were averaged. The noise contribution to our MSD data 
was estimated from the MSD of an immobile sample, resulting in Ay'^ld'l^ « 4 x 10“^, i.e. a factor of about 2.5 times 
smaller than the Ay^(<) values measured at short times. 

Simulations. In our molecular dynamics simulations, a 50:50 binary Yukawa fluid of large and small spheres with 
size ratio 1.2 is investigated. The model parameters have been reported earlier [26l Eli 07] . Our simulations have been 
performed for samples consisting of N = 12800 particles and having dimensions = 26.66ds, Ly = 13.33ds, = 

53.32dg. We work in the NVT ensemble using periodic boundary conditions, the temperature being controlled by a 
Lowe thermostat [Hj. The mode-coupling critical temperature of the system is Tc = 0.14. The system is equilibrated 
at T = 0.15 and then instantaneously quenched to T = 0.10, where it is aged for 10^ At this time, the walls 

are generated by freezing particles at 0 < z < 2ds and L^—2ds < z < [1H|. Stress is applied by pulling one wall 

at a constant force Fq in the x direction. For each applied stress, runs over 24 independent replicas of the system 
were averaged. Similar to the experiments, the dynamics were measured in a slice at the centre of the volume having 
thickness 13.3d5 and distance about 18ds to the walls on each side. 
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